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I. CURRENT PROGRESS 

During the first two months of the contract, the effort has proceeded 
at a low level while programs were being integrated and debugged. 

John Watts, SSL-MSFC has contributed greatly to program preparation. 

The transformations to appropriate coordinate systems have been 
checked. The Proton Spectrometer (PS) is rigidly mointed to the MDA 
structure with the cone opening in the -Y direction of the MDA coordinate 
system. The column vector representing this direction is (0, -1, 0). 

The SKYBET tape furnishes Euler angles to transform to the OWS 
coordinate system. This transformation accounts for bending at the 
flexible Airlock-OWS coupling. The OWS direction cosines of the PS 
ax' * are as follows: 

PSOWS1 = - cos El sin E3 cos E2 - sin El sin E2 

PSOWS2 ^ - cos El cos E3 

PSOWS3 = - cos El sin E3 sin E2 + sin El sin E2 


where 


(El, E2, E3) are the Euler angles to transform between the 
MDA and OWS systems. 


The PS axis must be related to earth-oriented coordinates of date 
for the pitch-angle distribution analysis and the east-west asymmetry 
analysis. A series of transformations are required. 



E • D • B • A • C • 


( PSOWS1 \ 
PSOWS2 ) 
PSOWS3 J 


Here (PS1, PS2, PS3) are the desired direction cosines of the PS axis in 
the rotating earth coordinate system at the time of the measurement. 
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Matrix C transforms from the OWS system to the ECI (Earth 
Centered Inertial, of date 1950.0) system. The Z-uxts of the ECI 
system points along the north pole of the spin vector. The X-axis is 
toward the vernal equinox. In terms of the SKYBET-furnished Euler 
angles, (Ul, U2, U3), The C matrix is as follows. Note that 
CU1 = cos Ul and SU1 = sin Ul, etc. 

/CU1 • CU2 CU1 • SU3 • CU2 + SU1 • SU2 
C - I - SU3 CU1 • CU3 

V CU3 • SU2 CU1 • SU3 • SU2 - SU1 • CU2 


- SU1 * SU3 • CU2 + CU1 • SU2 
SU1 • CU3 

SU1- SU3- SU2 ♦CU1 • CU2 


Matrix A transforms a column vector in ECI (epoch 1950) to ECI 
(epoch 1973) coordinates. It is a precession correction. Matrix A is 
furnished by the Dudley Observatory, Albany, New York. It was checked 
against Transformation 33, pp. 418-419, of Methods of Orbit Determination, 
Escobal, John Wiley & Cons, Inc., N»w York, 1965, and found identical 
to within two units in the seventh place after the decimal. 


A = 


.99998 

42912 

-.00514 

04016 

-.00223 

46372 

.00514 

04016 

.99998 

67880 

-.00000 

57437 

. 00223 

46372 

-.00000 

57433 

.99999 

75032 


Matrix B corrects for nutation (Wobble)from ECI (epoch 1973) 
to ECT (Earth Centered True, epoch 1973). This small correction is 
based on observational data and may be computed only after the fact. 
Matrix B is furnished by the Dudley Observatory. 


.99999 .99968 -.00007 37907 -.00003 19975 

B = [ .00007 .37903 .99999 .99972 -.00001 09138 

.00003 .19983 .00001 .09111 .99999 99994, 
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Matrix D rotates the ECT coordinate system around the Z-axis 
so that the X-axis is moved from the vernal equinox to the Greenwich 
mer^ian. 




e = 99.6909833 + 3 6000.7689 • TU + .00038708 • TU 2 


TU = elapsed time in Julian centuries, or 
TU = (2441682.5 - 2415020. )/36525 


A half day has been subtracted to change from Greenwich noon, 

Jan 0, 1950, to 0 hr, 0 min, Jan 0, 1973. 

Matrix E corrects from ECT (1973.0) to the time of measurement. 

( cos wt sin wt 

- sin wt cos wt 

0 0 

t = time since 0 hr., 0 min., Jan 0, 1973 
w= .26251 61453 radians/hr 

The six rotations shown above convert the PS axis from the local 
MDA system to the rotating earth system. 

The counts in proton channels P1-P8 are screened for validity. 
Total counts fewer than 10 distributed among the 8 channels during a 
1.2 second period are discarded. Counts greater than 50 in a single 
channel indicate a probable error and that period is discarded. The 
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counts are normalized to counts per second and dead time c .'rrections 
are applied. 

In this analysis, the measured values will be compared to values 
derived from the Vette AP7 proton environment model. An attempt will 
be made to compare spectral data. However, preliminary results point 
to a pain change on the Instrument. If corrections aio not possible, 
then the integral flux above 50 MeV will be compared. 

In keeping with these objectives, the analysis finds the AP7 
omnidirectional fluxes, ties a range of assumed pitch angle distributions 
to it, folds in instrument response, and computes the values which would 
be seen by an idea’i-ed PS immersed in the AP7 model environment. 
These calculations are made at each measurement point, i.e., every 1.2 
seconds along the trajectory where sensible flux is encountered. The 
measured values and calculated values are summed over time, then 
ratiod to provide spectral correction factors. 

At the same time, the data are tabulated according to location 
(B, L), pitch angle (a), energy, assumed pitch angle distribution (<r), and 
orientation in the earth-fixed coordinate system (® L y» * L y). the 

aid of this data, detailed corrections to the AP7 model environment are 
derived. "Best fit" energy -dependent pitch angle distributions are also 
obtained. Finally, some information is derived concerning the east- 
west asymmetry, though this effect is blurred at higher energies due 
to acceptance in the backward cone of the PS above 100 - 150 MeV. 

In order to accomplish the above objectives, the analysis 
continues from the point at which corrections were applied to the 
measured c’ata. 

The Vette AP7 flux parameters, a and b, are found by table 
look-up using values of B and L from the SKYBET tape. 
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AP7 FLUX ( > E) * a e" b E 

This form is differentiated, and the resultant is integrated over energy 
ranges corresponding to those of the eight proton channels. The angle 
a between the local magnetic field line, B, and the axis of the PS 
primary cone is found using the scalar product. The angle o is used 
in a table look-up to determine that fraction of the AP7 omnidirectional 
flux which would be observed by an ideal PS. The measured angular 
response of the instrument is triangular, peaked on the axis and going 
to zero at + 22. 5°. The assumed pitch angle distribution is Gaussian: 

-(90 - a) 2 

?V2 

e 

tr V 2 w- 

The idealized instrument reading is obtained by numerical integration 
and condensed into tables for this program. Note that a- determines 
the angular distribution of the flux. A set of 10 values are chosen for <r 
so that parallel calculations may be made. The best value will later 
be determined for each energy group, or the total flux. 

The flux values computed above will be denoted by $ . The 
measured fluxes will be denoted by «t>*. Each quantity will be further 
identified by two sets of subscripts. 

The first set of subscripts is described as follows. 

♦fjkl ’ ^ijklm 
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where 


i » B' (note prime) 
j * L' (note prime) 

k a a 
1 > E 

m a a 

Figure 1 shows that rectangular boxes in B - L coordinates would be 
wasteful of computer space, or would have large variations in flux 
level within a box. Therefore a transformation is applied to obtain 
B' - L' coordinates which are better suited here. 


/ - . 12942 
\- 1. 19972 

e a -3° 23’ 



2374 

9073 



The B' box limits are 0-. 004, .004-. 01, .01-. 02, and >.02. The L' 
box limits go from 0 to . 7 in steps of . 1, then >. 7. 

The box limits in the a dimension are 0 - 40, 40 - 60, 60 - 70, 
70 - 75, 75 - 80, 80 - 85, and 85 - 90 . Angles greater than 90° are 
reflected back to the first quadrant. 

The box limits in the energy dimension are identical with 
those of the PS. 

Thus, the measured fluxes, <t>*, are tabulated in a 4-dimensional 
table with 4x8x7x8 = 1792 subdivisions. The calculated fluxes, <J>, 
have an extra dimension, pitch angle distribution <r, and 10 parallel 
calculations are made for this quantity. 
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Flux units are 


p( E > 50 M cV I 

cm a - see 



Figure 1, Flux Contour Map of the Sky lab Orbit in B - L Space 
235 n. m. - 50° 

, p\G® ^ 



A second set of subscripts rind tables is necessary frr the east- 
west asymmetry study. These subscripts are: 

♦• tjl nl n2 ; nl n2 

where 

i ■ B' m * o’ 

J « L' nl * $ LV 

1 « E n2 » * LV 

The a subscript, k, is replaced by two subscripts, nl and n2, identifying 
equal solid angle boxes In the Local Vertical system. Here nl labels 
the polar angle w .h boundaries at 0, 60, 90, 120 and 180°, while n2 
labels the azimuthal angle with limits at 330, 30, 90, 150, 210, and 
270°. The solid angle Is subdivided into 4 x 6 = 24 equal size boxes. 
Here, the measured fluxes, <t>*, are tabulated in a 5-dimensional table 
with 4x8x8x4x6= 6144 subdivisions. 
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Values of $>*, <J>* , and <t> are determined for e.ich measurement 
(1.2 seconds) and are accumulated in separate tables according to the 
two sets of subscripts described abova. Auxiliary tables are compiled 
to show the number of entries in each box of each table. 

The tables described above are too lengthy to keep in fast 
storage. Therefore <!>*, the 10 <t>'s, and the two sets of subscripts will 
be placed in a buffer. When the buffer is full it is dumped into mass 
storage. After the input tape is processed, the values are brought back 
into fast core and processed. The subscripts, maximum value, and 
quantity represented are listed below: 
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E 

m 

10 

L’ 

n2 

6 + LV 

B’ 

nl 

49 lv 


k 7 


The measured and calculated data are sorted into tables as 
described above. The cumulative tables are stored on tape in the 
following order. 

Number of Entries Tables 

1. Number ol entries per bin for the pitch-angle tables 
(8 E's, 8 L's, 4 B's, 7a 's = 1792 words) 

2. Number of entries per bln for the east-west tables 
(8 E's, 8 L's, 4 B's, 6 azimuthal angles, 

4 polar angles - 6144 words) 

Pitch-Angle Tables 

3. Measured pitch-angle, $*, 1792 x 2 words (<& , <l>* 1 2 ) 

4-13. Calculated pitch-angle for 10 asrumed distributions, <X> , 

10 x 1792 x 2 words 
East-West Tables 

14. Measured east-west table, <t>*, 6144 x 2 words 

15-24. Calculated east-west tables for 10 assumed distributions, 

$ , 10 x 6144 x 2 words 
m 


The total table storage is 182528 words. 

Values are accumulated over each 12 hour period, then added to 
previous accumulations. Within each box, the mean values are: 




<*> = 7T £ $ » 


<t> 


1 V 4>* 2 

n £ - J 
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where n is the number of entries in that box. The standard deviation 
of ♦* is estimated from: 


(S.D.) 2 - i £(♦*- ± £(<t>* 2 - <*>* 2 ) 


Here the subscripts have been omitted for brevity. Note that summation 
over some of the subscripts is feasible before determining: means and 
standard deviations. 


It is planned that the following results will be computed 
following the processing of each 12 hour tape. 


A pitch angle distribution analysis will be made. The flux 
values will be summed over B' and L\ These sums are called 3* and 
S for the sums of <$>* and <t>, respectively. Then 


rt(E,a «-) 


SME, a ) 

a ,<r) 


F*(E)/g(0) t(0 d/3 

FtErftrFTrnmM 


Here, R(E, a , or ) is the ratio of observed to computed flux. The 
parameter <r labels the 10 assumed pitch angle distributions. F*(E) 
is the true omnidirectional flux in the energy band E. F(E) is the AP7 
model omnidirectional flux in the same band. The top (bottom) integral 
specifies the fraction of omnidirectional flux counted by the real 
(idealized) PS. The limits on the integrals involve and j but are 
complex because of the wide acceptance cone of the PS. The quantity 
t(/9) concerns the angular response of the instrument. The quantity of 
g(/3) is the true (but unknown) pitch angle distribution, while f(0, <r ) is 
one of the assumed distributions, 
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(90° -0) 2 

‘ 2 «r 2 

f(/9, ) = — ■ -- ■ 

The values to be extracted from the second equation above are 

F*(E) / F( E) and or ; R(E, a , a ) is known. Consider the following: 

the limits on the integrals are identical, as is tne PS angular response, 

t(/9). Therefore, if a choice of f(0) can be found that matches the true 

pitch angle distribution, g(/9), then the integrals will be identical and 

will cancel. In effect, the a dependence would vanish. A plot of 

R(E, c , <r ) versus a would be a horizontal line with ordinate F*(E)/F(E). 

For this reason, the best value of a, for a given energy, will be 

determined by doing at least squares fit to a horizontal line in 

R(E, a , <r ) versus a . The first one or two a values (0 to 40° or 60°) 

may have to be omitted due to sparse statistics. Alternatively, flux 

weighting may be employed. 

In the above manner, the best pitch angle distribution and ratio 
of true to AP7 omnidirectional fluxes will be obtained for each energy 
channel. The latter is simply the desired spectral correction factor 
to th . P7 model averaged over the Skylab orbit. 

The alternate tables using the second set of subscripts mentioned 
previously will be used to study the east-west effect. Ratios of the 
measured to calculated flux will be used here also because the orientation 
* *he spectrometer may not sample different dii actions equally in a 
.e time. Thus, the measured flux entering wr.h polar angle between 
• -j - 90° (nl = 2) and from the east 330 -30° (n2 = 1, divided by the 
corresponding calculated flux is: 

A(nl =2. n2 = 1) = 1) 
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Here the subscripts on B', L', E and <r are omitted. Particular values 
of these subscripts may be examined, or summation over desired ranges 
may be made prior to this calculation. A similar quantity may be made 
for flux in the opposite direction: 

A(nl - 3, n2 = 4) * gjj* . 

The east-west ratio for this choice is 


H 


_A(2 f 1) 

e-w ■ xry f "*T 


Fractional standard deviations may be computed for these values in the 
manner illustrated previously. 

The spectral correction factors, pitch angle distributions, and 
east-west factors will be computed over an increasing data base as the 
mission progresses. Increased accuracy should be reflected in decreased 
fractirnal standard deviations. 
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II. 


PROBLEM AREAS 


Spectral data may not be available due to gain changes in the PS. 
An attempt will be made to correct for gain changes. If this correction 
is not possible, the integral flux above 18 MeV will be analyzed. 
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ID. FUTURE WORK 

Over the next period, production runs will be made for the 
proton data. The codes will be modified to treat the electron and 
03 channel data. 
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